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We know very little about primordial curvature perturbations on scales smaller than about a 
Mpc. Measurements of the /x-type distortion of the CMB spectrum provide the unique opportunity 
to probe these scales over the unexplored range from 50 to 10 4 Mpc -1 . This is a very clean probe, 
in that it relies only on well-understood linear evolution. Also, just the information about the low 
multipoles (I ~ 100) of /i is necessary. We point out that correlations between /i-distortion and 
temperature anisotropies can be used to test Gaussianity at these very small scales. In particular 
the pT two-point correlation is proportional to the very squeezed limit of the primordial bispectrum 
and hence measures /jv£, while /i/i is proportional to the primordial trispectrum and measures tnl- 
We present a Fisher matrix forecast of the observational constraints on /]££ and stress that a cosmic 
variance limited experiment could in principle reach A/Jy£ ~ 0(1O~ 3 ). 



The initial conditions of our universe can be described 
in the simplest case just by specifying at some early 
time the probability distribution function of the adiabatic 
mode on a wide range of scales. Since then most of these 
scales have evolved in a complicated and often non-linear 
way until today. Therefore it is very hard to find accurate 
probes of these initial conditions, with two notable excep- 
tions: large scale structures today, which still have not 
entered a fully non-linear regime, and the anisotropies 
of the Cosmic Microwave Background (CMB) radiation. 
Because of complicated non-linear dynamics on the one 
side and Silk damping on the other, both probes are use- 
ful only at scales of order a Mpc or larger. Hence it is 
very interesting to find ways to investigate and constrain 
scales outside of the {10~ 4 — l}Mpc _1 window, which 
we have explored so far. 

In [T] it was shown that[T7], due to the dissipation of 
acoustic waves, the spectral distortion of the CMB can be 
used to constrain the (integrated) primordial power spec- 
trum in the approximate range 50 < fcMpc < 10 . This 
mechanism provides us with the probably unique chance 
to probe primordial perturbations that by now have been 
completely erased by Silk damping and swamped by com- 
plicated gravitational dynamics. Although we know that 
the primordial perturbations in the CMB/LSS window 
are close to Gaussian and hence well described by just 
the power spectrum, we know near to nothing about the 
statistical properties at much smaller scales. 

In this paper we show that, besides probing the in- 
tegrated power spectrum, the CMB spectral distortion 
is a potentially powerful tool to constrain deviations 
from Gaussianity at scales 50 < A: Mpc < 10 4 . Corre- 
lations between the /x-type distortion and temperature 
anisotropies of the CMB provide a direct measurement of 
the primordial bispectrum (in the squeezed limit), while 
self-correlations of /i-distortion measure the primordial 
trispectrum. It is important to stress that this probe is 
particularly clean in that it only relies on well understood 
CMB physics and is largely unaffected by gravitational 
non-linearities. 



The rest of the paper is organized as follows. We 
start by reviewing the relevant aspects of CMB physics. 
Then, following [THS], we derive an analytic estimate 
for the /^-distortion which we use to compute two-point 
self-correlations and cross correlation with temperature 
anisotropies. Finally we present a Fisher matrix fore- 
cast together with some consideration of observational 
prospects and we conclude with a summary. 

CMB DISTORTION 

At early times (z 3> 10 6 ) the universe is well described 
by a hot photon-baryon plasma. The number density of 
photons n(v) per frequency interval is given to very high 
accuracy by the black body spectrum n(y) = (e x — 1) 
where x = hv / '(kgT). The equation that describes the 
subsequent evolution of the photon number density is 
the full Boltzmann equation (when restricted to Comp- 
ton scattering this is known as the Kompaneet's equation 
[3]). This equation has three interesting regimes. Before 
z ~ Zfj,,i e 2 x 10 6 any energy released into the photon- 
baryon plasma is quickly thermalized by elastic and dou- 
ble Compton scattering (e _ 7H — > e~+2j), which are still 
very efficient [5]. The end result is again a black-body 
spectrum with now a higher T and a larger total number 
of photons N. After z^i double Compton scattering (as 
well as bremsstrahlung) becomes less efficient and the to- 
tal number of photons is approximately frozen [18j . For 
z /j,J ~ z ^ z \j,,ii w ith Zpj = 5 x 10 4 , equilibrium is still 
achieved after an energy injection due to elastic Comp- 
ton scattering, but since this process does not change 
the number of photons the end result is a Bose-Einstein 
distribution [BJ, i.e. n(v) — [e x+fJ "( x " > ~ l] where fj, is a 
frequency dependent chemical potential (rescaled by fc^T 
so that it is dimensionless) . The Kompaneets equation 
shows that /i(x) deviates from a constant only at very 
low frequencies, i.e. fj,(x) — noe~ Xc / x , with x c ~ 5 x 10 -3 . 
Henceforth we approximate \i as constant, which is valid 
everywhere except deep in the Rayleigh- Jeans tail (y — > 
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0). Finally for z < z^j even Compton scattering is not 
efficient enough to establish kinetic equilibrium between 
matter and radiation. The distortion created after this 
moment is known as y-type and is relevant e.g. for the 
Sunyaev-Zel'dovich effect [7J. Of course this is a simpli- 
fied picture since there is no sharp transition between one 
regime and the next. For the purpose of analytical esti- 
mates we will take the period responsible for the creation 
of (U-distortion to be z„ j < z < Za,% with the numerical 
values given above. As we will see, due to a logarithmic 
dependence on the size of this interval, changing these 
values by factors of order unity will not alter the main 
results. It should be clear though that for precise predic- 
tions one needs to study the system numerically. 

We will be interested in the energy injection coming 
from the dissipation of acoustic waves of the adiabatic 
mode (Silk damping) as these re-enter the horizon and 
start oscillating. Other sources of distortion are present 
(e.g. adiabatic cooling [2]) and the physics of the system 
is very rich. Our working assumption here is that either 
all other sources lead to a smaller and therefore negligi- 
ble distortion, as it is the case if the primordial power 
spectrum is not too red tilted, or that all other relevant 
effects are understood with a high enough precision to 
be subtracted off leaving the /^.-distortion caused by Silk 
damping as the only signal. 



Ai-DISTORTION 

In this section, following [IH5] , we derive a formula that 
relates the late time /x-distortion to the primordial power 
spectrum. Using the Bose-Einstein distribution plus the 
fact that the total number of photons is constant, for 
an amount of energy (density) released into the plasma 
SE one finds that p ~ 1ASE/E. Hence, let us estimate 
the energy injection due to damping of acoustic waves. 
The energy density of a density wave is given bv|19j Q = 
p(Sj(x) 2 ) p c 2 B I '(1 + Cg), with c s the sound speed, p the 
density and S the dimcnsionlcss amplitude of oscillations 
averaged over a period (indicated by () p to differentiate 
it from the quantum/ensemble average ()). Since at this 
time the universe is dominated by radiation we take p — 
p.y and c 2 ~ 1/3. Then one has 

SE F"' 1 d Q 1 
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We can use the transfer function (see e.g. [9]) 
A 7 (fc) ~ 3cos(fcr) e- fc2/fc ° , 
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Figure 1: The figure shows the power spectrum with Silk 
damping as function of logfc. The dotted, dashed and dot- 



dashed lines are A H e 



-2k A /k A D 



2 x 10 b 



5 x 



10 4 and zl = 1100 respectively. The red area on the right 
indicated by p is the difference of the power spectrum between 
and z M ,/. Once integrated over logfc this gives the p- 
distortion. For comparison on the left we have highlighted 
the scales probed by LSS and CMB anisotropics. 
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We then have 



(5 7 (a0 2 ) p = 1.45 
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x(A 7 (fc 1 )A 7 (A: 2 )) pe l ( £l+fc ~ 2 )- s , 

where 1Z describes curvature perturbations. Finally to 
account for the fact that p arises from a thermalization 
process, we use a top-hat filter in real space W(x), which 
smears the dissipated energy over a volume of radius 

ft s ~ DJ ' 

Summarizing, the deformation parameter p is related 
to primordial perturbations by 



p{x) 



4.6 
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x (cos (k\r) cos (k2r)) p 



where W(k) = 3fc -3 [sin(fc) — k cos(fc)] is the Fourier 
transform of the top- hat filter W(x) and k± = k\ ± 
k 2 . The quantum/ensemble average of p{x) gives the 
log-integral of the primordial power spectrum from 
k D (z^i) ~ 1.1 x 10 4 Mpc^ 1 to k D (zp,f) ^ 46Mpc~ 1 



(p{x)) ~ 2.3 / dlogfcA^(fc) 



(7) 
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where 



(K 2 ) = (2ir) 3 5 3 (K tot ) P n {k) , 
Pn{k) = 



27T 2 A|(fc) 



(8) 



To visualize this, in figure [l] we plot A 2 z (k)e^ 2k ' l l k o for 
kuiz^ i) (dotted line) and koiz^j) (dashed line) as well 
as their difference (red-filled region on the right), which 
quantifies the amount of dissipated energy and hence /z- 
distortion. 



TWO-POINT CORRELATION FUNCTIONS 

Let us now compute the two point correlation between 
a direction dependent /^-distortion /i(fi) and temperature 
anisotropy AT(h). To make contact with the way obser- 
vations are analyzed, we will decompose both signals in 
spherical harmonics, which are an orthonormal basis of 
functions on the sphere. We start with 



7 T = 
Hm — 
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where Ai(k) is the radiation transfer function. Through- 
out this paper we will use the Sachs- Wolfe approximation 
Ai(k) ~ jj(fcri)/3 with tl — 14Gpc the distance from 
last scattering and ji a spherical Bessel function. This is 
in principle valid only in the range 10 < I < 50, but for 
the purpose of the Fisher forecast it will be a reasonable 
estimate also for higher l's. Let us define in general 

((«L)*4w) = *«'U-cf , (io) 

then the TT correlation gives the well known result 



TT _ 2vr A n (k p ) 2 ^ 6.0 x 10~ 10 
1 ~ 25 1(1 + 1) ~ 1(1 + 1) 
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where k p = 0.002 Mpc -1 is a common pivot scale at 
which A^(fc p ) = 2.4 x 10~ 9 |8J. The /x-distortion can 
also be expanded as 

r d^k d^k - - 
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one can conveniently rewrite (12 1 as 



(7 3 k d^k 
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The Tfj, correlation is found to be 
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where we used the local bispectrum|20j 
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This formula can be simplified if one assumes a weak 
scale dependence. Taking A 2 z (k) — A 2 z (k p )(k / kp) 71 ' -1 
one finds 



kp 
2k, 



(n s - 1) log 
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1) close to the scale invariant 



so that b ~ 1 
limit n s — > 1. 

The fifi self-correlation has a Gaussian and a non- 
Gaussian contribution. For the former one finds 



l. Gauss 
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where in the last line we took An(kr>j) = Aji(k p ) and 
in order to get an upper bound k s = fc^j, which is 
the shortest scale that contributes sizably to dissipation. 
This scaling can be understood as follows. For a Gaussian 
field, consider the power spectrum (i.e. /^-distortion) at 
small scales kjjj. Its fluctuations in regions much more 
distant than fc^ ^ are completely independent. Hence, 
the expansion of these power fluctuations in multiples is 
independent of I, i.e. just white noise, for low multiples 
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<C kr>,f - Now let us consider the effect of smearing. 
Measuring the low /^-distortion multiples, we are effec- 
tively averaging over many small-scale (again of order 
knj) independent realizations, and hence we are sup- 
pressing the individual variance by N~ 1 / 2 . The number 
of kr> f modes in a spherical shell around the last scatter- 

2\ 



ing surface of thickness k s is given by N ~ k 3 D j/(k s r L 
Using the local trispectrum 

(K 4 ) = (2irj 3 S 3 ^K tot ^T NL 

x Pn{ki)Pn{k 2 )Pn{\ki + h\) + 11 perm's 



for the non-Gaussian contribution one finds 
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where b' is defined by 
V 



1(1 + 1) 



din k + dln k-d\n k% 



A R (fc_/2) 2 A R (fc + ) 2 A R (fc 3 ) 2 
A n (k P ) e 

lk 3 /k D 

i 



X 




I 








f 





(21) 



such that b' ~ 1 for a scale invariant power spec- 
trum. Notice that, for a given amount of primordial non- 
Gaussianity, the ratio of the non-Gaussian to Gaussian 
contributions to the two-point self correlation is much 
larger for fi than it is for T. In other words, Cf* is more 
sensitive in relative terms to non-Gaussianity than Cf T . 

To discuss the practical observability of the above sig- 
nals we need to consider experimental noise. For T this 
is negligible since we will only consider low multiples 
(I < 0(100)). To model the noise for fi we assume a 
Gaussian beam and use [11] 
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where is the sensitivity to /x and Z max is related to the 
beam size of the experiment. For example, for an exper- 
iment like PIXIE [TS] the beam size is 9fwhm = l°-6 

— 1/2 

which leads to l ma x — 84 and the sensitivity is ~ 
v4~7T x 10~ 8 . So a figure of merit to keep in mind is 



1.3 x 10~ 15 e' 2/ '" 



FISHER MATRIX ANALYSIS 

We are now going to use the above results to perform 
a Fisher forecast for the bounds that /^-distortion can 



put on /jvl local. The fact that we invoked the Sachs- 
Wolfe approximation instead of using numerical transfer 
functions should not alter much the estimate of the signal 
to noise ratio because the acoustic peaks in T mostly scale 
away. The lxl Fisher matrix 
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gives us the signal-to-noise ratio via S/N = JnlVF- For 
mild non-Gaussianity, Jnl J; 10 5 , we can estimate the 
noise in the measurement of each I as 



of 
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where we used that fj, and T instrumental noises are un- 

- as seen by 



correlated, Cj T > Cf 



and 
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comparing (19) with (23). Then the signal-to- noise ratio 



can be written as 
S 



(25) 



where (S/N)^ is the signal-to-noise ratio for the \x- 
distortion averaged over the whole sky. Apart from the 
numerical coefficient, this relation is to be expected. The 
long-scale fluctuations of \i around its average that corre- 
late with temperature fluctuations are generated by the 
non-Gaussian part of 1Z and hence are suppressed with 
respect to the // monopole by JnlTZ- Notice that reduc- 
ing the beam size of the experiment improves the signal 
to noise ratio only logarithmically. Plugging in numbers 
one finds the figure of merit 
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where A/i = 10 8 is the estimated one sigma error on 
the jti-distortion monopole of an experiment like PIXIE 



and b, defined in (18), is of order one if the primordial 
power spectrum is approximately scale invariant. 

There are at least two classes of models in which the 
fiT correlation could provide the strongest constraints on 
non-Gaussianity already with PIXIE's sensitivity. First, 
models in which the power spectrum grows at small scales 
|13) . since then b ^> 1. In most of these models b can be as 
large as 0(1O 2 ) leading to A/^£ ~ 0(10) for PIXIE. Sec- 
ond, models in which the bispectrum diverges faster than 
(23) the local template in the squeezed limit, i.e. (1Z 3 ) cx k^ a 
for &3 — > with a > 3 (see e.g. [TJ] for a phenomenolog- 
ical model). 

Before concluding, let us remark that while instrumen- 
tal noise can be improved, there is a lower bound for the 
noise imposed by nature, often referred to as cosmic vari- 
ance. In this respect the fiT cross correlations possess an 
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important advantage as compared with the temperature 
bispectrum, which is another probe of the primordial lo- 
cal bispectrum [21]. Cosmic variance, and more generally 
sampling variance, scales with the number of modes N as 
^y-i/2_ p or a p r j morc ii a i bispectrum of the local type it is 
useful to distinguish between long and short modes. The 
latter are the same for both observables, provided one 
measures the low multiples of fj, (I ~ 100). On the other 
hand, the number of short modes in the temperature bis- 
pectrum is about Z^ ax and due to Silk damping Z max can 
never be larger than a few thousand. In fj,T cross corre- 
lations the number of short modes is much larger and, as 
we said, can be estimated as k? D f /(k s rl 2 ) > 10 12 . This 
means that, for comic variance limited experiments, the 
bispectrum's sensitivity is A/jvl ~ C(5), while using fiT 
correlation one can in principle reach A/nl ~ O(10 -3 ). 



that CMB spectral distortion will one day tell us another 
little bit about the very early stages of our universe. 
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CONCLUSIONS 

Measurements of the /i-distortion of the CMB spec- 
trum offer the unique opportunity to probe primordial 
perturbations in the otherwise inaccessible range 50 < 
fcMpc < 10 4 . We have shown that the two-point self 
correlations of //-distortion anisotropies and cross cor- 
relations with temperature anisotropies provide a direct 
measurement of the primordial tri- and bispectrum in 
the squeezed limit, respectively. We have performed a 
Fisher matrix forecast for the bounds on /jvl local and 
found that a bound of |A/jvl| < 10 3 is achievable already 
with current technology [15] . It is worth stressing again 
that in the case of mild non-Gaussianity the monopole 
/i-distortion is expected to be larger than the anisotropic 



part (see (25) for a quantitative statement). 

It would be interesting to further study the perfor- 
mances of a dedicated experiment. We should stress that 
one should be careful in comparing this bound with those 
from CMB/LSS, since it applies to very small scales, 
around a kpc, which are completely unconstrained so far. 

We have also noticed that the /zT cross correlation is 
contaminated by a much smaller cosmic variance than 
the temperature bispectrum |22 and for an ideal exper- 
iment it could probe an /jvx much smaller than unity. 
In this respect it is important to stress that for all prac- 
tical purposes single-field inflation predicts a primordial 
bispectrum that vanishes in the squeezed limit. The well 
known consistency relation [16] /nl = n a — 1 is derived 
using comoving momenta. When comparing with obser- 
vations one should use physical momenta instead. Doing 
this leads to an additional term that exactly cancels the 
n a — 1. As long as one considers other sources of primor- 
dial local non-Gaussianity so that /jvl 3> n a — 1, this 
slow-roll correction can be safely neglected as we have 
done in the present derivation. 

Finally, given the pace with which observations have 
improved in recent years, there is good reason to hope 
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Strictly speaking this is not true deep in the Rayleigh- 
Jeans tail (y — > 0), but it is a good approximation for the 
range of frequencies relevant for observations. As we will 
see shortly this is quantify by the frequency dependence 
of the chemical potential. 

This formula contains the correct relativistic factor 
c s/(l + c s) an d f° r small c\ reduces to the non-relativistic 
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result used e.g. in [T]. For c 2 s — 1/3 The difference is ex- 
actly the factor 3/4 found in [12j . 
[20] To simplify our notation we omit the label loc in /jy£. 
Since local non-Gaussianity is the only form of non- 
Gaussianity we consider in this work, this should not 
lead to any confusion. 



[21] Remember though that the two observables constrain 
different scales and are therefore complementary rather 
than competing. 

[22] It would be interesting to check whether any astrophysi- 
cal foreground can contaminate the signal. 



